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We present preliminary measurements of branching fractions for the semileptonic decays B — » 
Dt~v t and B — ► D*t~V t , which are potentially sensitive to non-Standard Model amplitudes. 
The data sample comprises 232 x 10 6 T(4S) — > BB decays collected with the BABAR detector 
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at the PEP-II e+e" storage ring. We obtain B(B~ -» D°t~V t ) = (0.63 ± 0.38 ± 0.10 ± 0.06)%, 
B{B~ -* L L>*°r-77 T ) = (2.35±0.49±0.22±0.18)%, B(B° -» D + t~V t ) = (1.03±0.35±0.14±0.10)%, 
and B(B° -» D* + t~V t ) = (1.15 ± 0.33 ± 0.04 ± 0.04)%, where the uncertainties are statistical, 
systematic, and normalization, respectively. By combining B~ and B° results, we also obtain 
the branching fractions B{B -> Dt~V t ) = (0.90 ± 0.26 ± 0.11 ± 0.06)% and B(B -> D*t'V t ) = 
(1.81 ± 0.33 ± 0.11 ± 0.06)% (quoted for the B~ lifetime), with significances of 3.5cr and 6.2<r. 

Submitted to the 2001 Europhysics Conference on High Energy Physics, Manchester, England. 
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Semileptonic decays of B mesons to the r lepton — 
the heaviest of the three charged leptons — provide a 
new source of information on Standard Model (SM) pro- 
cesses 

dual 

, as well as a new window on physics beyond 
the SM [J, la, la, [3, l8| ■ In the SM, semileptonic decays oc- 
cur at tree level and are mediated by the W~ boson, but 
the large mass of the r lepton provides sensitivity to ad- 
ditional amplitudes, such as those mediated by a charged 
Higgs boson, H~ . Experimentally, b — ► ct~v t decays are 
challenging to study because the final state contains not 
just one, but two or three neutrinos. 

Branching fractions for semileptonic B decays to r 
leptons are predicted to be smaller than those for I = 
e, fi 9], but are still substantial compared to most 
hadronic B decays. A recent SM-based calculation Q 
predicts B(B° -> D+t-V t ) = (0.69 ± 0.04)% and 
B(B° D*+t-V t ) = (1.41 ±0.07)%; an inclusive calcu- 
lation gives B{B -> X c t-v t ) = (2.3 ± 0.25)%, where 
X c represents all final states resulting from the b — *■ c 
transition. Calculations [4J, |5|, |g, |7J, |8[ in supersymmetric 
models show that substantial departures from the SM de- 
cay rate could occur for B(B — > Dt~V t ), but that those 
for B(B — > D*t~v t ) are expected to be smaller. The 
interference with the SM amplitude can be constructive 
or destructive, depending on the value of (tan /3)/m,H , 
where tan (3 is the ratio of the vacuum expectation val- 
ues for the two Higgs doublets and tuh is the H~ mass. 

Theoretical predictions for semileptonic decays to ex- 
clusive final states require knowledge of the form factors, 
which parametrize the hadronic current as a function of 
q' 2 = (pb —Pd(*)) 2 - F° r light leptons (e, fi), there is ef- 
fectively one form factor for B — > D£~17i, while there are 
three for B — > D*l~vt. If a r lepton is produced instead, 
one additional form factor enters in each mode. The form 
factors for B — ► D^*'l~Vi decays involving the light lep- 
tons have been measured and have been discussed 
extensively in the theoretical literature. Heavy-quark- 
symmctry (HQS) relations [ll| allow one to express the 
two additional form factors for B — ► Z)(*V _ i7 T in terms 
of the form factors measurable from decays with the light 
leptons. With sufficient data, one could probe the addi- 
tional form factors and test the HQS relations. 

The first measurements of semileptonic 6-hadron de- 
cays to r leptons were performed by the LEP experi- 
ments [12j operating at the Z Q resonance, yielding an 
average [13] branching fraction Z?(&had — * Xt~v t ) = 



(2.48 ± 0.26)%, where &had represents the mixture of b- 
hadrons produced in Z° — + bb decay. 

We determine branching fractions of four exclusive 
decay modes Q: B~ -> D°t~V t , B~ -> D* t~V t , 
B° — > D + t~V t , and B Q — * D* + t~V t , each of which is 
measured relative to the corresponding e and \x modes. 
To reconstruct the r, we use the decays r~ — * e~V e v T 
and t~ — > pTV^VT, which are experimentally most ac- 
cessible. The main challenge of the measurement is to 
separate B — > D^*'t~V t decays, which have three neu- 
trinos, from B — * D^(~U£ decays, which have the same 
observable final-state particles but only one neutrino. 

We analyze data collected with the BaBAR detec- 
tor [ijl], at the PEP-II e + e~ storage ring at the Stanford 
Linear Accelerator Center. The data sample used in the 
analysis comprises 208.9 fb _1 recorded on the T(4S) res- 
onance, yielding 232 x 10 6 BB decays. The measurement 
uses all of the major detector subsystems: a charged- 
particle tracking system consisting of a 5-layer silicon ver- 
tex tracker and a 40-layer He-gas drift chamber (DCH); 
a quartz-bar Cherenkov particle-identification system; a 
CsI(Tl) crystal electromagnetic calorimeter for electron 
identification and photon energy measurement; a 1.5 T 
superconducting magnet; and a muon identification sys- 
tem in the magnet flux return. 

The analysis strategy is to reconstruct the decays of 
both B mesons in the Y(4S) — > BB event, providing 
powerful constraints on unobserved particles. One B 
meson, denoted Btag, is fully reconstructed in a purely 
hadronic decay chain. The remaining charged tracks and 
photons are required to be consistent with the products 
of a b — > c semileptonic B decay: a hadronic system, 
either a D or D* meson, and a lepton (e or fj,), either 
primary or from r~ — * t~viy T . Using the known to- 
tal four-momentum of the e + e _ collision, we calculate 
Pmiss = b(e + e~) -ptag -Pn(«) —pe] recoiling against the 
observed B ta , g + D^*'£ system. A large peak at zero in 
TO miss = Pmiss corresponds to semileptonic decays with 
one neutrino, whereas signal events form a broad tail out 
to m^ iss ~ 8 (GeV/c 2 ) 2 . To separate signal and back- 
ground events, we perform a fit to the joint distribution 
of w^gg and the lepton momentum (|p||), in the rest 
frame of the B meson. In a signal event, the observed 
lepton is the daughter of the r and typically has a soft 
spectrum; for most background events, this lepton typi- 
cally has higher momentum. 
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We reconstruct B tag candidates [l6[ in 1114 final states 
Stag — ► D^Y ± . Tag-side candidates are recon- 

structed in 21 decay chains, and the system can con- 
sist of up to six light hadrons (n ± , it , K ± , or K®). B tag 
candidates are identified using two kinematic variables, 
toes = \/s/4:- |pta g | 2 and AE = E tag - where 
^/s is the total e + e~ energy; |ptag| is the magnitude of 
the -Btag momentum; and -Bt ag is the Btag energy, all 
defined in the e + e~ center-of-mass frame. We require 
toes > 5.27 GeV/c 2 and \AE\ < 72 MeV, corresponding 
to 4<t (standard deviations). We reconstruct B tag candi- 
dates with an efficiency of approximately 0.3% to 0.5%. 

For the B meson decaying semileptonically, we re- 
construct D^*' candidates in the modes D° — ► A'~7r + , 

K-TT+TT , K-TT+1T+7T-, K°7T + 7T~; D+ -> A'-7T + 7r+, 

K-ir+ir+n , K° s ir+, K-R+-K+; D*° -> D°ir°, L>°7; 
and D*+ -> L> tt + , D + tt°. D (£>*) candidates are se- 
lected within 4cr of the D mass (D* — D mass differ- 
ence), with cr typically 5-10 MeV/c 2 (1-2 MeV/c 2 ). To 
ensure well-measured momenta, identified electron and 
muon tracks are required to have at least 12 hits in 
the drift chamber and not to be near the acceptance 
edges. Electron candidates must have lab-frame momen- 
tum |p e | > 0.3GeV/c; muon candidates must have an 
appropriate signature in the muon detector system, effec- 
tively requiring IpJ ^ 0.6 GeV/c. The energy of electron 
candidates is corrected for bremsstrahlung energy loss if 
photons are found close to the electron direction. 

We require that there be no charged tracks not asso- 
ciated with the -Btag, D^*\ or £ candidates. We com- 
pute -Bextra, the sum of the energies of all photon candi- 
dates not associated with the -B tag , D^*\ or £ candidates, 
and we require -B extra < 150-300 MeV, depending on the 
£>M channel. We suppress hadronic events and combi- 
natoric backgrounds by requiring |p m i SS | > 200 MeV/c 
and q 2 > 4 (GeV/c 2 ) 2 . If multiple candidates pass this 
selection, we select the candidate with the lowest value 
of B extra . To improve the to 2 ^^ resolution, we perform 
a kinematic fit to the event, constraining particle masses 
to known values and requiring tracks from B, D, and K® 
mesons to originate from appropriate common vertices. 

Figure [T] shows the distributions of for the four 

D^*'i channels, along with the projections of the max- 
imum likelihood fit to be discussed below. The large 
peaks at ~ are mainly due to B —* D^*H~Vi, 

which serve as normalization modes. The structure of 
this background is shown in the inset figures, which 
expand the region -0.4 < m 2 ^ < 1.4 (GeV/c 2 ) 2 . 
B — > D*l~T>t background is the dominant feature in 
the two D*£ channels (Figs. QJ,, c); the two Dl chan- 
nels (Figs. [Tp, d) are dominated by B — > D£~V( de- 
cays but also include substantial feed-down contributions 
from true D* mesons where the low-momentum 7r° or 
photon from D* — ► Dir° /j is not reconstructed. This 
feed-down is clearly visible for B — > D*£~Vt background, 
but affects B — > D*t~v t signal similarly, and both feed- 



down components (as well as smaller feed-up contribu- 
tions from B — > D(£~ /t~)T> into the D*£ channels) are 
included in the fit. Other sources of background include 
B — > D**(£~ /t~)V events (here D** represents charm 
resonances heavier than the Z?*(2010), as well as non- 
resonant D^mr systems); charge-crossfeed (which oc- 
curs when a B — > D(*'l~Vt event is reconstructed with 
the wrong charge for the B tag and meson, typically 
because a low-momentum is swapped between the 
-Btag and the D^); and combinatoric background. This 
last background is dominated by hadronic B decays such 
as B — > D^Di*\ in which one of the charm mesons pro- 
duces a secondary lepton, including r leptons from D s 
decay. 

To constrain background from B — > D**(£~ /t~)v de- 
cays, we use four control samples (one for each signal 
channel) in which an extra ir° meson is observed. Most 
of the D** background in the signal channels occurs when 
the 7T° from D** -► £)(*) 7T° is not reconstructed, so these 
control samples provide a good normalization of the back- 
ground source. -D** decays in which a ^ is lost do not 
have the correct charge correlation between the -B tag and 
, and decays with two missing charged pions are rare. 
The feed-down probabilities for the D**(£~ /t~)V back- 
ground are determined from simulation, with uncertain- 
ties in the D** content treated as a systematic error. 
However, the control samples reduce our sensitivity to 
the details of this model. 

We perform a relative measurement, extracting both 
signal B — > D(*>t~V t and normalization B — > D^£~T7( 
yields from the fit to obtain the four branching ratios 
R(D°), R(D+), R{D*°), and R(D*+) where, for exam- 
ple, R(D*°) ee B(B- -» D*°t-V t )/B(B- -> D*°£-p e ). 
Here, £ represents only one of e or /i; however, both 
light lepton species contribute statistically to the denom- 
inator. Signal and background yields are extracted us- 
ing an extended, unbinned maximum likelihood fit to 
the joint (m 2 ^, |p||) distribution. The 18-parameter 
fit is performed simultaneously in the four signal chan- 
nels and the four D** control samples. In each of 
the four signal channels, we describe the data as the 
sum of seven components (shown in Fig. [I]): Dt~T> Ti 
D*t-V t , D£~ut, D*£~V tl D**(£-/t-)V, charge cross- 
feed, and combinatoric background. The four D** con- 
trol samples are described as the sum of five compo- 
nents: D**(£~ /t~)V, D£~V<t, D*£~Vt, charge crossfeed, 
and combinatoric background. Probability distribution 
functions (PDFs) are primarily determined from simu- 
lated event samples; however, the parameters describing 
the dominant feed-down component — D* feed-down into 
the D£ channels — are determined directly by the fit. 

We perform two fits, one in which all four signal yields 
are allowed to float independently, and a second, B~~B° 
combined fit, in which we constrain 

and R(D* + ) = R(D*°). The fit results are summarized 
in Table HI and the to 2 ^^ projections of the constrained 
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FIG. 1: Distributions of events and fit projections in m^ iss for 
the four signal channels: D*°£,_D°£, D* + £, and D + £. (The 
fit shown incorporates the B~-B° constraints.) The normal- 
ization region m^iss ~ is shown as an inset in each figure. 
The fit components are: D*t~V t (green), Dt~V t (purple), 
D*rv t (blue), Drv e (red), D**{i~/T~)V (yellow), com- 
binatoric (grey, below dashed line), charge crossfeed (grey, 
above dashed line). 



parametrization (2% to 12%); the composition of com- 
binatoric backgrounds (2% to 11%); the mixture of 
D** states in B -> D**£~V e decays (0.4% to 6%); the 
B -> D*e~V e form factors (0.1% to 1.8%); and the tt° 
efficiency, which affects the D* — ► D feed-down rate 
(0.4% to 1.0%). Uncertainties on the w^iss resolution 
for B — > D*t~vi events and on the B — > D£~V( form 
factors each contribute less than 1%. The net systematic 
uncertainty on i? associated with fit yields is given by 
(Ai?/i?)fit for each channel in Table [U Uncertainties on 
R propagated from the ratio of efficiencies for signal and 
normalization modes are typically small due to cancella- 
tions, and include the limited statistics in the simulation 
(1.1% to 1.5%) and systematic errors related to detector 
performance. The latter are determined by studying the 
efficiency of track and neutral reconstruction and parti- 
cle identification performance in control samples in data 
and contribute less than 0.2% each, except for and ^ 
identification, which contribute 0.5% to 0.7% each, and 
are larger because the lepton momentum spectrum differs 
between the signal and normalization processes. Finally, 
the uncertainty on B(t~ — > l~Viv T ) 13J contributes 0.2% 
to all modes. The net systematic uncertainty on R due 
to the efficiencies is given by (AR/R) e in Table Q] 

These results are preliminary. We estimate that uncer- 
tainties in R due to modeling of bremsstrahlung radiation 
are at or below the 1% level, but they have not been ex- 
plicitly included in the results. While B — > D^*'£~Vt 
decays are modeled with HQET-based form factors [l9| 
including recent ex per imental measurements [l(| , we cur- 
rently use ISGW2 [20j to model signal decays. 

We determine the statistical significance of the sig- 
nals from ^/2A(ln£), where A(ln£) is the change in 
log-likelihood between the nominal fit and the no-signal 
hypothesis. The total significance is determined in a sim- 
ilar manner, by modifying the likelihood function to take 
into account systematic uncertainties from the fit. Ta- 
ble U gives both significances for each channel. 

We have presented preliminary measurements of the 
decays B — ► Dt~V t and B — ► D*t~17 t , relative to the 
corresponding decays involving light leptons. We obtain 
R(B -► Dt-V t ) = (40.7 ± 12.0 ± 4.9)% and R(B -> 
D*t-V t ) = (31.0 ± 5.7 ± 1.8)%, where the first error is 
statistical and the second is systematic. Normalizing to 
known branching fractions [13j, we obtain 



B(B^Dt-V t ) = (0.90 ±0.26 ±0.11 ±0.06)% 
B(B — > D*t V t ) = (1.81 ±0.33 ±0.11 ±0.06)%, 



fit are shown in Fig. [TJ 

Systematic uncertainties on R associated with the fit 
are determined by running ensembles of fits in which in- 
put parameters are distributed according to our knowl- 
edge of the underlying source, and include the PDF 



where the third error is the uncertainty on the normaliza- 
tion branching fraction, and where results are expressed 
for the B~ lifetime. The significances of the signals 
are 3.5cr and 6.2tr, respectively. The measurement of 
B — > D*t~V t is consistent with a preliminary Belle mea- 
surement (2lj : the measurement of B — > Dt~V t is the 



7 



TABLE I: Preliminary results from fits to data and associated uncertainties: the columns are the signal yield (JV s j g ), the yield 
of normalization B — > D^t~Vi events (7V n orm), the ratio of signal and normalization mode efficiencies (esig/eWm), the relative 
systematic error due to the fit yields ((AR/R){[ t ), the relative systematic error due to the efficiency ratios ((AR/R) e ), the 
branching fraction relative to the normalization mode (7?), the absolute branching fraction (B), and the total and statistical 
signal significances (a to t and <7 sta t). The first two errors on 7? and B are statistical and systematic, respectively; the third error 
on B represents the uncertainty on the normalization mode The last two rows show the results of the fit with the B~-B° 
constraint applied, where B is expressed for the B~ . 



Mode 


N ■ 


■^V norm 


^■sig / ^norm 


(AR/R) Rt 


(AR/R) E 


R B a tot 










[%] 


[%] 


[%] [%] (<T B tat) 


B~ -> 


D u t~V t 33.1±19.6 


346.7±23.0 


1.85 


15.2 


1.5 


29.5±17.4±4.5 0.63±0.38±0.10±0.06 1.7 (1.7) 


B- -> 


D*°t-V t 95.9±19.8 


1628.6±63.5 


0.98 


9.4 


1.4 


36.2± 7.5±3.4 2.35±0.49±0.22±0.18 5.3 (5.8) 


B° -► 


D+t~V t 23.0±7.9 


149.9±13.3 


1.83 


13.4 


1.7 


48.6±16.7±6.6 1.03±0.35±0.14±0.10 3.3 (3.5) 


B° -» 


D* + t~V t 16.2±7.3 


481.8±25.5 


0.93 


3.4 


1.5 


21. 4± 9.7±0.8 1.15±0.52±0.04±0.04 2.7 (2.7) 


B -» 


Dt~V t 64.9±19.1 


496.3±26.4 


1.85 


12.0 


1.3 


40.7±12.0±4.9 0.90±0.26±0.11±0.06 3.5 (3.8) 


B -> 


D*t~V t 105.3±19.4 


2109.4±68.0 


0.93 


5.7 


1.2 


31.0± 5.7±1.8 1.81±0.33±0.11±0.06 6.2 (6.5) 



first evidence for this mode. These results are about la 
higher than predictions based on the Standard Model, 
but, given the size of the uncertainty, there is still room 
for a non-SM contribution. 
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